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' Abstract. Neutron stars, with their strong surface gravity, have interestingly short timescales for the 

^ I sedimentation of heavy elements. Recent observations of unstable thermonuclear burning (observed 

O ■ as X-ray bursts) on the surfaces of slowly accreting neutron stars (< 0.01 of the Eddington rate) 

' motivate us to examine how sedimentation of CNO isotopes affects the ignition of these bursts. We 

further estimate the burst development using a simple one-zone model with a full reaction network. 
We report a region of mass accretion rates for weak H flashes. Such flashes can lead to a large 
reservoir of He, the unstable burning of which may explain some observed long bursts (duration 
- 1000 s). 
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^ INTRODUCTION 

O _ 

An ionized plasma in a gravitational field develops an electric field sufficient to levitate 
Ph| the ions and ensure overall charge neutrality. When there is more than one species of ion 

P ■ present, the ions will experience a differential force: lighter ions float upward (defined 

by the local gravitational field) and heavier ions sink downward. 
^ ! Accreting neutron stars, with their strong surface gravity ^ 2.0 x 10^^ cms^^, are an 

ideal place to look for the effects of sedimentation. The sedimentation of heavy elements 
and resulting nucleosynthesis in the envelope of isolated neutron stars cooling from 
birth was first described by Rosen [1,2] and has been studied in detail by Chang and 
I Bildsten [3, 4]. For accreting neutron stars, the rapid stratification removes heavy nuclei 

from the photosphere for accretion rates M < lO^^^Moyr^^ [5]. Deeper in the neutron 
star envelope, the differentiation of the isotopes can alter the nuclear burning, namely 
unstable H/He burning and rapid proton-capture process (rp-process), that powers Type 
I X-ray bursts. 

Recent long-term monitoring of the galactic center with BeppoSAX led to the discov- 
ery of nine "burst-only sources" [see 6, and references therein]. These sources did not 
have persistent fluxes detectable with the BeppoSAXfW¥C, thus must have extremely 
low accretion luminosities (Lx < 10^^ ergs s^^ ~ O.OlLEdd)- If the accretion is not con- 
centrated onto a small surface area, so that the local accretion rate is m ^ M/ (AnR^), 
then the sedimentation timescale for CNO nuclei, defined as the time required to move 
a scale height relative to the center of mass of a fluid element, is less than the accretion 
timescale, defined as the time for a fluid element to reach a given depth. An estimate of 
this accretion rate was noted earlier by Wallace et al. [7]. As a result, sedimentation in 



X 



the accreted envelope must be considered in treating the unstable ignition of hydrogen 
and helium for these low-m sources. At somewhat higher accretion rates, Lx ~ O.OIL^^, 
several bursts have been observed with longer durations, < 1000 s, that are intermediate 
between mixed H/He bursts and superbursts. 

Motivated by these observations, we explore the unstable ignition of hydrogen and 
helium at low accretion rates and pay particular attention to the regime where the mass 
accretion rate is less than the critical rate needed for stable H burning [8]. We also study 
how the sedimentation of CNO nuclei affects the unstable ignition of H and He in an 
accreted neutron star envelope. For a full reporting of these results, please see [9]. 



THE ROLE OF SEDIMENTATION ON IGNITION CONDITIONS 

To study the accumulation of H and He, we construct models of a plane-parallel atmo- 
sphere is in hydrostatic equilibrium, and for simplicity we neglect the terms coupling 

thermal and particle diffusion [see 10, and references therein]. We follow the treatment 
of Burgers [11] and construct for each species an equation of continuity and momentum 
conservation. Collisions are described by a resistance coefficient Kgi — nsnt{OstVst), with 
Vst the center-of-mass relative velocity between particles of types s and t. 

As an illustration of how sedimentation changes the structure of the envelope, consider 
a trace species, labeled 2, in a background of a species labeled 1 , i.e., ^ ^2 and species 

1 has velocity wi = 0. The electric field is eE — m-^gA\j{Zi + 1) where the electrons are 
non-degenerate and eE = m^^A\IZ\ where the electrons are degenerate. Substituting E 
into the equation of motion for species 2 then determines the sedimentation velocity, 
w'sed = wi= «2(^2Wug — ZieE^jKyi. We chose the sign of Wged to be positive if species 

2 moves downward. 

Using the Stokes-Einstein relation to determine K\i from the drag coefficient for 
a liquid sphere ([12]), and a nonrelativistic electron equation of state, we find the 
sedimentation velocity of a trace nucleus [see 13] 

3 r0-3 (A2Zi-Z2Ai)A0i _i 
Wsed = 2xl0 gu-^ ^3-03 cms , (1) 

where we use the common shorthand = g/llO^^ cms'^), = r/(10^K), and 
P5 = p/(10^gcm~^). In a pure H plasma, the sedimentation velocity^ is greater than the 
mean velocity m = m/p for m < 400 g cm~^ s~ ^ T^-^p^-^ [A^ —Zq)/!!^^; under conditions 
at which H ignites (Ty ~ 5; ps ~ 4) this corresponds to m < 0.02mEdd and m < O.OSwEdd 
for a trace ^He and ^^C nucleus, respectively. Here mgdd = 2mpc/[(l + Xyi)(7thR] — 
8.8 X 10^ g cm^ s^ is the local Eddington mass accretion rate for a solar composition, 
<7th is the Thomson scattering cross section, and Xh is the hydrogen mass fraction. 

We numerically solve for the evolution of the accreting envelope prior to unstable 
ignition by coupling the diffusion equations to the equations for temperature and heat 



Throughout this paper we use a Newtonian metric and assume a neutron star mass 1.4 Mq and radius 
10 km, so that ^14= 1.9. 
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where dy — —pdz is the column density, Cp is the specific heat at constant pressure and 
Vad = {d\nT/d\nP)s. For detailed description of the numerical method in solving these 
equations, please refer to [9]. 

We take the composition of matter accreted from the companion star to be roughly 
solar composition in and ^He, and distribute the remaining mass evenly between ^^C 
and 1^0: X(1h) =0.71, X(4He) =0.27, X(i2c) = 0.01, and X(160) =0.01. Below the 
fuel layer are ashes from previous X-ray bursts; we set the ash composition to ^^Zn, 
consistent with the findings of recent one-dimensional calculations of repeated X-ray 
bursts [14, 15]. We use a radiative-zero outer boundary, and for the inner boundary, at 
the bottom of the ash layer, we set the flux to a constant value, Fc = 0.1 MeV m/mu [16]. 
This flux is set by reactions in the deep crust — electron captures, neutron emissions, and 
pycnonuclear reactions. 

We explore the evolution of the accreting neutron star envelope for accretion rates 
10~^-0.6mEdd (10-5 X 10"^ g cm^ s^ for accretion of a solar composition mixture). For 
each accretion rate, we ran calculations both with and without isotopic sedimentation, 
and we follow the thermal and chemical evolution of the mixture until the envelope 
becomes thermally unstable. We determine this ignition point using a simple one-zone 
criterion, {derMc\/dT)p > {dEcooi/ dT)p, where Cnuci is the nuclear heating rate and 
£cooi = pKT/y^ is an approximation of the local cooling rate. The mass fraction of 
H and He at the column where either H or He unstably ignites is plotted in fig. 1, as a 
function of mass accretion rates. When sedimentation is included, the abundance of H 
at the base of the accreted envelope is depressed. It is tantalizing that the accretion rate 
at which mixed H/He ignition occurs is increased by a factor of 2 when sedimentation is 
taken into account, and we speculate that this might alleviate the discrepancy between 
the predicted transition in burst duration [8] and recent observations [see, for example 



When the temperature in the neutron star envelope is sufficiently low, the CNO cycle be- 
comes temperature dependent; as a result, the ignition of H becomes thermally unstable 
at low accretion rates [8]. As a first investigation of the unstable ignition of H when the 
atmosphere is stratified, we perform several calculations of the burst nucleosynthesis. 
We approximate the cooling by a one-zone finite differencing over the envelope. 



where £cooi = pKTy^ and we evaluate Cnuci from a reaction network. Included in this 
network are 686 isotopes covering proton-rich nuclei up to Xe [see 18]. The reaction 
rates are taken from the compilation reaclib [see 19, 20, and references therein], and 
consist of experimental rates and Hauser-Feshbach calculations with the code NON- 
SMOKER [20]. The initial temperature, column density, and composition are taken from 
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FIGURE 1. Mass fraction of hydrogen at the column where either H (open symbols, solid lines) or He 
(filled symbols, dotted lines) unstably ignites, as a function of mass accretion rates. We show results for 
which sedimentation is ignored (circles), and for which it is included (squares). From [9] and reproduced 
by permission of the AAS. 

the values of the quasi-static calculation at the ignition point. For the unstable ignition of 

H, the temperature in the accreted envelope is set by the flux emergent from the deeper 
ocean and crust. As a result, the accretion rate at which the burst behavior changes will 
depend on assumptions about the heating in the crust. 

We shall first describe the outcome of these one-zone calculations (eq. [3]). The left 
panel of Figure 2 shows the post-ignition evolution for accretion rates 9.1 x 10^ rh^dA 
(80 g cm^^ s^^; solid lines), 1.1 x 10^-^ m^AA (100 g cm^^ s^^; dotted lines), and 
2.3 X 10^^ WEdd (200 g cm^^ s^^; dashed lines). We plot three different quantities: the 
evolution of temperature {top panel), the heat flux, Fcooi = y£cool> normalized to the 
accretion flux (middle panel), and the mass fractions of hydrogen and helium (bottom 
panel). We then repeat this calculation at higher mass accretion rates (Fig. 2, right panel), 
5.7 X 10~^ WEdd (500 gcm^^ s^^; solid lines) and 1.1 x lO^^mEdd (10^ gcm^^ s^^; dot- 
ted lines). 

It is immediately clear that there are two very different outcomes for the H ignition. 
At lower accretion rates, the rise in temperature following H ignition is sufficient to 
trigger a vigorous He flash with a decay timescale ~ 10 s (Fig. 2, middle-left panel). At 
higher accretion rates, however, the flash is too weak to ignite He. As an interpretation 
of these results, recall that the ignition curve for the 3a reaction has a turning point 
at jturn ~ 5 X 10^ gcm^^, which terminates the unstable branch [21, 22]. At the higher 
accretion rates (Fig. 2, right panel) the ignition of H occurs at j < jtum> and the local 
rise in temperature does not trigger unstable He ignition. In fact, the temperature rise is 
not even sufficient to initiate convection, as the radiative temperature gradient needed to 
carry Fcooi is dXnT /dXnP ~ 0.25 < {dXwT / d\nP)s. In contrast, at lower accretion rates 
(Fig. 2, left panel), the ignition of H occurs at _y > jtum, and the rise in temperature will 
ignite the triple- a reaction. 

If subsequent H flashes do not ignite the underlying He, then a large layer of nearly 
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FIGURE 2. One-zone burst calculation following unstable H ignition. Left: for three mass accretion 

rates, ih/ihEdd = 9.1 x 10^^ (solid lines), 1.1 x 10^' {dotted lines) and 2.3 x 10^^ (dashed lines), re- 
spectively. Right: for two mass accretion rates, m/mEdd = 5.7 x 10^^ and 0.011. Top panel: temperature 
evolution; Middle panel: the ratio of one-zone cooUng flux to the accretion flux; Bottom panel: mass frac- 
tion of hydrogen (thin lines) and ^He (thick lines), respectively. Diffusion and sedimentation is included. 
From [9], and reproduced by permission of the AAS. 



pure He will accumulate. Because our system is not burning H steadily, the temperature 
in the He layer is colder than if the burning were in steady-state. As a result, a large 
He layer should accumulate. To get the thermal structure of an accumulated pure He 
layer, we integrated equations (2), with d jdt 0, from the column where H ignition 
occurred, and set the temperature there to that found in the accumulating model. For 
w/mEdd = 5.7 X 10^ and 0.011, the He ignition column density varies from y 
2 X 10^-3 X 10^^ gcm~^ (recurrence time ~ 0.06 — 19 yr) for Fcmu/rh varies from 1.0 
to 0.1. Our one-zone approximation gives a cooling time T 400 s[};/(10^'^ gcm~^)], 
which roughly agrees with the long burst duration. 



CONCLUSION 

Using a simplified numerical model of an accreting neutron star envelope that allows for 
differential isotopic velocities, we have undertaken a first study of the effects of sedimen- 
tation and diffusion on the unstable ignition of hydrogen and helium on the surfaces of 
accreting neutron stars, and have investigated the outcome of unstable hydrogen ignition 
using simple one-zone models. Our principal conclusions are: 

1. The effect of sedimentation changes the conditions at H/He ignition even for 
m>0.1mEdd- 

2. There is a range of accretion rates for which unstable H ignition does not trigger 



a He flash. This range depends on the flux from the deep crustal heating and the 
degree of settling of He and CNO nuclei. We speculate that successive weak H 
flashes can lead to the accumulation of a large reservoir of He; this may explain the 
long bursts observed from some sources. 
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